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ABSTRACT 

We present results from a study of the kinematic structure of star-forming galaxies 
at redshift z ~ 3 selected in the VVDS, using integral- field spectroscopy of rest- 
frame optical nebular emission lines, in combination with rest-frame UV spectroscopy, 
ground-based optical/near-IR and Spitzer photometry. We also constrain the under- 
lying stellar populations to address the evolutionary status of these galaxies. We infer 
the kinematic properties of four galaxies: VVDS-20298666, VVDS-020297772, VVDS- 
20463884 and VVDS-20335183 with redshifts z = 3.2917, 3.2878, 3.2776, and 3.7062, 
respectively. While VVDS-20463884 presents an irregular velocity field with a peak 
in the local velocity dispersion of the galaxy shifted from the centre of the galaxy, 
VVDS-20298666 has a well-resolved gradient in velocity over a distance of ~4.5 kpc 
with a peak-to-peak amplitude of v =91 km s _1 . We discovered that the nearby galaxy, 
VVDS-020297772 (which shows traces of AGN activity), is in fact a companion at a 
similar redshift with a projected separated of 12 kpc. In contrast, the velocity field of 
VVDS-020335183 seems more consistent with a merger on a rotating disk. However, 
all of the objects have a high local velocity dispersion (a ~ 60-70 km s _1 ), which gives 
v/a<l. It is unlikely that these galaxies are dynamically cold rotating disk of ionized 
gas. 

Key words: galaxies: evolution - galaxies: formation - galaxies: kinematics and 
dynamics - instrumentation: adaptive optics - galaxies: high redshift - galaxies: star- 
burst 



1 INTRODUCTION 

To understand the assembly of stars in galaxies, it is crucial 
to characterize the star formation rate and the total stellar 
mass already formed as a function of the dynamical mass 
(the mass of the halo including dark matter). Determining 
this relation at intermediate and high redshifts allows us to 
explore whether we are seeing the initial build-up of stellar 
mass in a fully- formed dark halo, or if the stellar mass to 
total mass ratio is fairly constant compared to lower-redshift 
samples (as might be the case if much of the stellar mass 
formed very early, or dark haloes were assembled through 
dry mergers triggering little further star formation). The 
redshifts range z ~ 2 — 4 is therefore of great interest as 
massive galaxies are thought to undergo their most active 
period of star formation (e.g, Hopkins & Beacom 2006). 
In the past, deep multi- wavelength surveys (e.g., 
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GOODS - IGiavalisco etTaE] [20041 Hubble De ep Field - 
|Williams et al.|1996| DEEP2 - [Faber et al.|2007| have high- 
lighted the large diversity in morphology and stellar pop- 
ulations among high redshift galaxies. Many studies using 
UV-IR photometry and long-slit spectroscopy, mainly fo- 
cusing on the Lyman-break galaxy (LBG) populations (e.g. 
|Erb et al.|2003 2004 2006a| and complemented by various 



studies using lensed objects (e.g., Lemoine-Busserolle et al. 



2003 Bunker et al. 2000), have helped to understand the 
global characteristics of z ~ 2 — 4 galaxies. However, with 
the limited information available from long-slit spectroscopy 
it is difficult to answer keys questions on the dynamical na- 
ture of high redshift galaxies. It is plausible that some may 
be perturbed by bursts of star formation and merging to lie 
outside the low-redshift morphological classification scheme 
(the familiar Hubble sequence) . To properly explore whether 
this is the case, or instead if most star- forming galaxies have 
stable disk kinematics, requires studying the resolved kine- 
matic structure through "3D" spectroscopy with Integral 
Field Units (IFUs). The recent advent of high resolution 
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Table 1. SINFONI Observations 



Galaxy 


VVDS-ID 


RA (J2000) 


DEC (J2000) 


z(l) 


Exp. Time (2) 


Seeing (3) 


Airmass (4) 


Dates of Obs. (U.T.) 


VVDS-3884 


020463884 


02 27 01.23 


-04 21 04.4 


3.2780 


3.5 Hrs 


0'/38 


1.28 


07 Sept. 2005 


VVDS-8666 


020298666 


02 26 18.32 


-04 20 32.7 


3.2805 


5 Hrs 


/ /54 


1.12 


06,08 Sept. 2005 


VVDS-7772 


020297772 


02 26 18.23 


-04 20 34.1 


_t 


5 Hrs 


/ /54 


1.12 


06,08 Sept. 2005 


VVDS-5183 


020335183 


02 25 33.71 


-04 15 41.4 


3.6993 


6 Hrs 


0'/38 


1.07 


13,15 Nov. 2006 



t This nearby galaxy to VVDS-8666 was discovered in our IFU data cube to share the same redshift. The original photometric redshift 
from VVDS/CFHT-LS was z p h ot = 1.747 with broad minimum in x 2 . It was not observed spectroscopically with VIMOS as part of the 

VVDS survey. (1) redshift estimated from rest-frame UV spectrum obtained with VIMOS. (2) Total observing time usable for each 
target (900s for each individual exposures). (3) median seeing estimated from the PSF stars. (4) average airmass of all the exposures. 



IFUs on 8m-class telescopes has enabled the study of high- 
redshift galaxies on scales of a few kiloparsecs in seeing- 



limited conditions (e.g. |Bunker et al.| 2QQ4; Sm ith et al.|2 QQ4) 
and even down to sub-kiloparsec scales using adaptive op- 
tics and the magnification due to gravitational lensing (e.g., 
Stark et al. ^008). Very recent studies using Integral Field 



Spectroscopy have started to give an insight into the kine- 
matics properties and physical phenomena at play in high 
redshift galaxies, such as merging and galactic winds associ- 



ated with strong star formation (|Genzel et al.|2006; 


Law et 


al. 2007a, 2009; Wright et al. 2007, 2009; Forster Schreiber 


et al. 2006; 


Bouche et al.| 20071 Genzel et al. I2008| |Nes- 


vadba et al. 


|2008||van Starkenburg et al.|2008|). There is 



some controversy from this early work about the properties 
of the dynamical structure of these distant galaxies (particu- 
larly the proportion of rotationally supported gaseous disks 
versus objects kinematically dominated by high velocity dis- 
persions). However, there is a consensus from the work so 
far that at z > 2 galaxies typically exhibit higher veloc- 
ity dispersions (relative to ordered rotation) compared with 
lower-redshift samples. 

In this paper, we present results on the kinematic prop- 
erties of a sample of four galaxies at 3 < z < 4 selected in the 
VVDS Deep data (I < 24; |Le Fevre et al.|2005] ), using IFU 



i^-band spectroscopy with VLT/ SINFONI. These spectra 
are sensitive to the rest-frame optical, in particular the neb- 
ular emission lines from gas ionized by hot stars in regions 
of recent star formation. Hence we can determine the cur- 
rent star formation rate, and the spatial distribution of star 
formation, as well as measuring the kinematics of the gas. 
We are able to study the stellar populations through CFHT 
photometry from both the VVDS and CFHTLS surveys, and 
also imaging with the Spitzer Space Telescope using IRAC 
and MIPS, complemented by UV rest-frame spectroscopy 
from VLT/ VIMOS. The results presented here are part of a 
study to construct samples representative of the global inter- 
mediate and high-z population by selecting galaxies only on 
the basis of their magnitude (Iab < 24.75). Two companion 
papers, Lemoine-Busserolle et al. (2009b, in revision) and 
Queyrel et al. (2009, in revision) are devoted to the study of 
intermediate redshift galaxies (1 < z < 1.5) observed with 
VLT/SINFONI during the same observing runs. 

In Section[2]we describe our sample, observational strat- 
egy and the data reductions techniques. In Section [3] we ad- 
dress the nature of the stellar population of the galaxies, 
using broad-band photometry to constrain properties such 
as stellar mass, age and star formation rate. In Section [4] we 
investigate the kinematic structure and dynamical proper- 
ties inferred from the integral field spectroscopy. Finally in 



Section [5] we summarize our results and explore the nature 
of these LBG galaxies at z = 3 — 4. We compare our findings 
with previous studies of kinematic properties of z ~ 2 — 3 
galaxies. 

We assume a cosmology with Q = 0.3,A = 0.7 and H 
= 70 km s _1 Mpc -1 throughout, and all magnitudes are on 
the AB system ( |Oke fe Gunn|1983 ). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample Selection 

Our IFU observations targetted a sample of 3 galaxies se- 
lected in the 02h "deep" (17.5 < Iab < 24.0) field of the 



2005) 



VIMOS VLT Deep Survey (VVDS; [Le Fevre et al 
The z > 3 objects presented here were selected on the basis 
of redshifts derived from the VLT/ VIMOS rest-UV spec- 
tra, such that the expected wavelengths of the rest-optical 
H^/[OIII] emission lines would be clear of bright OH sky 
lines in the near-IR. 



2.2 SINFONI Observational Strategy 

The near-IR spectroscopic observations of the sample were 
acquired with the 3D-spectrogra ph SINFONI (SINg l e Fain t 
Object Near-IR Investigation 



Eisenhauer et al. (2003)) 



at ESO-VLT during two 4-nights runs, on September 5-8, 
2005 (ESO run 75.A-0318) and on November 12-15, 2006 
(ESO run 78.A-0177). We used the largest-field mode of 
SINFONI, to enable good sky subtraction without the need 
to do offset sky observations. This mode has a spaxel scale 
of / /l25 x Of! 25 pixels (with coarser sampling in y), lead- 
ing to a field-of-view of 8 /; x 8". Our observations used the 
i^-grism, spanning 1.94 — 2.46 /xm with 2.45 A pixels, and 
are detailed in Table [I] Unresolved sky lines had a mea- 
sured width of ~ 7 A (FWHM), indicating a spectral reso- 
lution Rs = A/AAfwhm ~ 3300 (or 90km/s velocity reso- 
lution) . Conditions were photometric and the median seeing 
for each objects is in Table ^ Each target was acquired 
through a blind offset from a nearby bright star. Individ- 
ual integrations were 900 s, and we read out the array in 
multiple non-destructive read mode to reduce the readout 
noise to ~ 7 e~ and become background-limited. Between 
readouts the telescope was nodded by w 4", positioning the 
galaxy in opposite corners within the 8" x 8" SINFONI 
field-of-view. Two sets of observations located around the 
upper corner (A) and the lower corner (B) were obtained 
per object. This observational procedure allows background 
subtraction by using frames contiguous in time, but with the 
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Figure 1. Two-dimensional spectrum of the nebular emission lines of VVDS-8666 and VVDS-7772. The x-axis is wavelength (increasing 
to the right) and spans the range A = 20800 — 21620 A. The y-axis is spatial distance, with a length of 6'ZO displayed. 



galaxy in different locations. Moreover, the target was never 
located exactly at the same position on the detector, a min- 
imal sub-dithering of 0.3" was required in order to minimize 
instrumental artifacts when the individual observations are 
aligned and combined together. The target was reacquired 
from the offset star every hour. The blind offset stars were 
also used to monitor the PSF. The total on-source integra- 
tion times are listed in Table [T] SINFONI was used in its 
seeing-limited mode for two galaxies (VVDS-020298666 in 
September 2005, and VVDS-020335183 in November 2006 
- hereafter abbreviated to VVDS-8666 and VVDS-5183), 
and in closed-loop NGS-AO (Natural Guide Star - Adap- 
tive Optics) mode for the observations of VVDS-020463884 
(hereafter called VVDS-3884) on 07 September 2005, which 
had a suitably bright K = 10.7 nearby star 27" 4 away. The 
Strehl ratio achieved ait 2.2 Lim was 15%. 



2.3 SINFONI Data Reduction and Flux 
Calibration 

Data reduction has been performed with the ESO-SINFONI 



pipeline (version 1.7.1, Modigliani et al.|[2007| ), IRAF and 
custom IDL scripts. The individual 900s exposures were run 
through the pipeline, which flat-fields the frames using an in- 
ternal lamp, corrects for hot pixels using dark frames, wave- 
length calibrates each individual observation from the spec- 
tra of reference arc lamps, and then reconstructs and spa- 
tially registers the three-dimensional datacubes (with the 
sky background still present) for each exposure after apply- 
ing a distortion correction. Our observing strategy was im- 
plemented to maximize on-source integration and therefore 
no offset blank-sky frame was obtained. As the field of view 
used with SINFONI is larger than the spatial extent of the 
galaxies, we were able to obtain a clean sky background sub- 
traction by stepping the galaxies between two well-separated 
points within the field of view, using the approach and IDL 



script presented in Davies (2007). Sky background subtrac- 



tion was done directly on the cube, using the individual 
cubes at position B as a sky for the individual cube at po- 
sition A and vice versa. After sky subtraction, individual 
cubes were aligned in the spatial direction by relying on the 
telescope offsets from the nearby bright star (used as a refer- 
ence for the blind acquisition) and then combined together, 
rejecting remaining cosmic ray strikes to produce the final 
reduced cube for each object. 

A flux calibration is required in order to derived abso- 
lute parameters (e.g. star formation rate) from the uncali- 
brated flux measured in emission lines. Each science obser- 
vation has been accompanied by the observation of a telluric 



standard star at similar airmass of B spectral type and mag- 
nitude K = 6 — 7. Integration times were 2 — 6 s for these 
standard star observations, and the cubes were reduced in a 
similar way to our science exposures. We extracted ID spec- 
tra of the stars (summing all the flux within an aperture 
of diameter 5 resolution elements) which were used both to 
perform the flux calibration of our galaxy cubes, and to cor- 
rect them from telluric absorption lines. We assume that the 
standard stars are well fitted by a pure blackbody curve (a 
good approximation for these B stars). Starting from the 
blackbody temperature of the standard star and its magni- 
tude in if -band, we calibrated the sensitivity function using 
IRAF. The galaxy cubes were then divided by the calibra- 
tion curve appropriate for that observation to produce a flux 
calibrated spectrum, corrected for atmospheric absorption. 
At the centre of the if -band, the conversion was approxi- 
mately 1 count/s = 2 x 10 -17 erg cm -2 s _1 . 

We measured the noise in each data cube by taking a 
slice in wavelength close to the position of the line emission. 
We then determined the RMS counts in blank areas uncon- 
taminated by emission from the target galaxy in the region 
where all the dithered pointings overlapped. In computing 
the flux uncertainty in an emission line from the extracted 
spectrum, we multiplied the measured pixel-to-pixel noise 
in the data cube by ^/7rr^ pat d\, where r spa t is the radius in 
pixels of the extraction aperture, and d\ is the wavelength 
extent (again in pixels) over which the line flux is measured. 
For the galaxies VVDS-5183, VVDS-8666 and VVDS-3884 
we used spatial apertures of r spa t = 8 pixels (r = 1//0), 
and for the compact companion VVDS-02029772 (hereafter 
abreviated to VVDS-7772) we used a smaller aperture of 
Tspat — 4 pixels (r = (Z/5). In all cases, we used an extrac- 
tion of width 8 pixels in wavelength (20 A or 300km/s), 
which was ^ FWHM of the lines. The typical 1 a noise was 
0.2 — 0.3 x 10 _17 ergcm -2 s _1 for our extraction aperture. 

The blind offset stars were reduced in the same way 
than the telluric standards. We produced collapsed 2D im- 
ages using the SINFONI pipeline to allow to the extent of 
the PSF to be estimated from the offset stars. 



3 SPECTROPHOTOMETRY PROPERTIES 

3.1 VVDS-7772: A companion galaxy of 
VVDS-8666 at z = 3.29 

Our IFU spectroscopy of VVDS-8666 reveals a companion 
galaxy with [O III] emission at a similar redshift. The [O III] 
emission of the companion is blue-shifted by 307 ± 34km/s 
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Table 2. Nebular emission lines measurements from SINFONI ID spectra 



Galaxy 


A(H/3) (1) 


z H( 3 (2) 


AA(H/3) (3) 


f(H/3) (4) 


A([Om] A5007) (1) 


z [0 in] ( 2 ) 


AAQOiii] A5007) (3) 


f([Om] A5007) (4) 


/([Om] A5007) 




/A 




FWHM / A 


/10" 17 erg cm- 2 s- 1 


/A 


FWHM / A 


/10" 17 erg cm- 2 s- 1 


//(H/3) 


VVDS-3884 


20797.81 


3.2782 


21.52 


5.67 ±0.28 


21417.50 


3.2776 


22.21 


16.17 ±0.34 


2.85 


VVDS-8666 


20866.36 


3.2923 


18.43 


4.44 ± 0.20 


21487.68 


3.2917 


14.96 


7.10 ±0.20 


1.59 


VVDS-7772 








< 0.3 (2<r) 


21468.79 


3.2878 


11.27 


2.58 ±0.14 


> 9 (2 a) 


VVDS-5183 


22876.63 


3.7058 


16.60 


7.40 ±0.30 


23563.85 


3.7063 


18.82 


24.10 ±0.50 


3.26 



(1) Wavelength of nebular emission line. (2) redshift estimated from rest-frame optical nebular emission line. (3) FWHM of emission 
line. (4) Emission line flux with uncertainties. Non-detections are recorded as 2 a upper limits. 
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Figure 2. One-dimensional UV-rest frame spectra of VVDS- 
3884,VVDS-5183 and VVDS-8666. Vertical dotted lines indicate 
the wavelengths of the strongest spectral features 



relative to the targetted galaxy (see Figs. [T] & |3j. The 
companion is included in the VVDS photometric catalog 
as object VVDS-020297772 (hereafter VVDS-7772), with 
Iab = 25.0, and is undetected in [/-band. The photometric 
redshift is only z p hot = 1.7478, very different from our spec- 
troscopic redshift of z spec — 3.2878, which may be due to an 
unusual SED (perhaps contaminated by the presence of an 
AGN, see Section [33). 



The spatial separation of the galaxies is l'/6 =b 0'/2 from 
the line emissions in the 3D IFU cube, which corresponds to 
a projected distance of 12kpc at z — 3.29 (the nominal sep- 
aration in the broad-band images in worse seeing is l'/9). If 
we assume that these two galaxies are gravitationally bound, 
and the velocity offset is due to orbital motion and not out- 
flows etc., the enclosed mass is given by v 2 r/G, and we can 
set a lower bound on the orbital radius, r, by using the pro- 
jected separation (12kpc) and also the resolved component 
of the orbital velocity along the line of sight (307 km s -1 ). 
This yields M syste m > 2.6 x 10 11 M©. 



3.2 Rest-Frame UV Spectroscopy 

Optical spectra have been obtained with the VIMOS in- 
strument installed at ESO/VLT as part of the VVDS survey 
(Fig. [2|. The spectral coverage is 5600 A < A < 9250 A which 
corresponds approximately to rest-frame UV wavelengths 
1300 A < A < 2150 A for the z w 3.3 galaxies, and 1190 A < 
A < 1970 A for VVDS-5183 at z « 3.7 (which encompasses 
Lyman-al216A in this high redshift case, but not in the 
other galaxies). The resolution is R s ~ 230. Of the three 
galaxies, VVDS-8666 has a relatively featureless spectrum 
with only O 1 1302 A/ C 11 1304 A absorption prominent. 
VVDS-3884 has strong absorption features from the high- 
ionization interstellar medium (ISM) lines C iv 1548,1550 A 
and Si iv 1394,1504 A (with rest-frame equivalent widths 
of Wo = 4.5 A and 3.0 respectively), as well as weaker 
C 11 1334 A and Si 11 1526 A absorption. The high-ionization 
lines are blueshifted with respect to the [O 111] redshift by 
« lOOkms -1 , with the low ionization lines blueshifted by 
« 300kms _1 . There is a hint of He 11 1640 A emission in 
the ID spectrum. In contrast, VVDS-5183 has more promi- 
nent low-ionization ISM absorption lines than the high- 
ionization C iv & Si iv. All the ISM lines are blueshifted 
by ~ 500kms _1 relative to the [O 111] redshift. This galaxy 
exhibits weak Lyman-a emission of 6 x 10 -18 ergs cm -2 s _1 
that is redshifted by 200kms _1 relative to [O 111], with a 
rest-frame equivalent width of Wo « — 7 A. 



3.3 Rest-Frame Optical Spectroscopy and 

Metallicity Estimates from Nebular Lines 

Integrated rest-frame optical ID spectra were extracted from 
the SINFONI data cube for each object. From these we 
measured line fluxes and velocity widths of the integrated 
line emission (see Fig. [3). We can use the [O iii]/H/3 ra- 
tio to investigate any possible contribution of an AGN to 
the nebular spectrum. We follow the classification proposed 
by |Lamareille et ah] |2009 ) . This classification is based on 
the distribution of galaxies of known type in the 2dFGRS 
survey, with respect to their [O m]/H/3 ratio. The classifi- 
cation as a star-forming galaxy is secure for log([0 iii]/H/3) 
lower than 0.4, and the classification as an AGN is secure for 
log([0 iii]/H/3) greater than 0.7. It is impossible to classify 
with certainty galaxies lying in the region between these two 
limits, but we know that this region is dominated by star- 
forming galaxies (« 60% for log([0 m]/H/3) lower than 0.6). 
With a log([0 m]/H/3) of 0.46 and 0.51, respectively, VVDS- 
3884 and VVDS-5183 are candidate star-forming galaxies. 
The nebular spectra of these two galaxies are probably pro- 
duced in hot H 11 star-forming regions. However, there is also 
a chance (~ 40%) that they are instead dominated by an 
AGN (i.e. Seyfert 2 galaxies). With log([0 m]/H/3) of 0.20, 
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Figure 3. SINFONI integrated ID spectra of the nebular emis- 
sion lines, R(3 to [O in] AA4959,5007(indicated by vertical dashed 
lines). At the bottom of each panel, a scaled version of the sky 
spectrum is shown (dotted line). Galaxy VVDS-5183 is shown 
in the top plot, with VVDS-8666 in the middle panel (with the 
spectrum of the nearby companion, VVDS-7772, shown offset by 
10 -17 above); the bottom panel shows the spectrum of VVDS- 
3884 (smoothed with a boxcar of 15 pixels; the unsmoothed data 
is shown offset above). 




10 

log(M*/M ) 

Figure 4. This is the stellar-mass-metallicity relation for 
the galaxies VVDS-3884 (diamond), VVDS-8666 (triangle) and 
VVDS-5183 (square). The solid curve shows the relation at z ~ 
( Tre monti et al.|2004|, the long-dashed curve shows the relation 
at z ~ 2 (|Erb et al.||2006a[), th e short-dashed curve the relation 
at z ~ 3 iMaiolino et al.|2008| ). 



VVDS-8666 is without doubt a star- forming galaxy. Finally, 
the Hf3 line from the companion (VVDS-7772) is undetected, 
unlike in VVDS-8666, and we derive an upper-limit for H(3 
emission by collapsing the data cube around the wavelength 
inferred from companion redshift from [O in] , and perform- 
ing aperture photometry and background subtraction on the 
2D line map. We quote a 2 a upper limit from the local back- 
ground noise in the 3D data cube. We can thus compute a 
lower limit for log([0 iii]/H/3) which is 0.93. With such a 
large [O m]/H/3 ratio, VVDS-7772 is probably dominated 
by an AGN. 

Gas-phase oxygen abundances may be also estimated 
from the [O m]/H/3 ratio (see Table [2} using the relation 



provided by Liang et al. (2006), which was calibrated on 
Tremonti et al.| fl2004p metallicities in the SDSS DR4 sam- 
ple. We find 12+log(0/H) of 8.57 ± 0.02, 8.66 ± 0.02, and 
8.56 ± 0.02 for VVDS-3884, VVDS-8666, and VVDS-5183, 
respectively. 

Figure [4] shows the position of these galaxies in the stel- 
lar mass-met allicity plane. In this figure, we also show the 

2 ( |Erb et al.||2006a| ), 



( |Tremonti et al.||2004|), z 



and z rsj 3 (Maiolino et al. 2008) relations, the last two 
being rescaled to Tremonti et al.| (2004) metallicities using 
equations from Kewley &; Ellison] ( 2008 ) . They fall between 
the z ~ 2 and z ~ 3 relations. We note however that the 
[O m] /Hp ( |Liang et al.||2006 ) calibration relies on the un- 
derlying relation between gas-phase oxygen abundance and 
ionization degree in the parent sample. Thus, this calibration 
might overestimate metallicities for galaxies at high redshift, 
provided that they are likely to show higher ionization states 
and higher [O iii]/H/3 ratios for lower metallicities than the 
ones observed in the local universe. Compared to the local 
mass-metallicity relation, the galaxies of our sample show a 
—0.47 dex mean systematic shift in met allicity. 
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Table 3. Photometry 



Galaxy 


Bab (1) 


Vab (1) 


Rab (1) 


Iab (1) 


uab (2) 


9AB (2) 


TAB (2) 


iAB (2) 


zab (2) 


Jab (3) 


Kab (3) 


3.6 Mm (AB) (4) 


4.5 /xm (AB) (4) 


VVDS-3884 


24.65±0.09 


23.77±0.06 


23.83±0.06 


23.39±0.08 












24.04±0.31 


23.29±0.20 


22. 0±0.lt 


22.8±0.3 


VVDS-8666 


25.6±0.19 


24.68±0.10 


24.03±0.09 


23.67±0.12 


27.06±0.34 


25.32±0.06 


24.26±0.03 


23.82±0.02 


23.72±0.06 




22.83±0.18 


22.1±0.1 


21.9±0.1 


VVDS-7772 


27.53±0.82 


26.30±0.34 


26.06±0.41 


25.01±0.32 


> 27 


26.31±0.19 


25.74±0.13 


25.52±0.13 


> 27 




22.57zb0.20* 






VVDS-5183 


> 27 


25.37±0.24 


23.95±0.07 


23.68±0.10 




26.17±0.22 


24.15±0.04 


23.77±0.03 


23.45±0.08 






22.7±0.2 


22.8±0.3 



t Confused source; * The K-band magnitude is anomalously bright for VVDS-7772 
The columns are as follows: (1) CFHT Legacy Survey observations (2) TERAPIX CFHTLS-T0003 observations (3) SOFI/NTT 

observations (4) IRAC photometry 



Table 4. Stellar population properties 



Galaxy 


M v 


Mi 


M* (a) 


Age (b) 


Luv (c) 


SFR UV (d) 


SFRsed (e) 


E(B - V) st ar (f) 


VVDS-3884 


-22.40 


-23.78 


1.29 


718 


23.85 ±1.38 


20 ± 1.1 


20.2 





VVDS-8666 


-22.86 


-23.68 


1.63 


29 


19.68 ±1.60 


16 ± 1.3 


77.9 


0.35 


VVDS-7772 


-23.12t 




2.85t 


32lt 


3.09 ± 1.17 


2.6 ±0.97 


ot 


ot 


VVDS-5183 




-23.19 


1.51 


4.79 


30.63 ±0.85 


25 ±0.7 


532 


0.40 



t The if-band magnitude, from which My is derived, is anomalously bright for VVDS-7772. We use a post-starburst SED with no dust to 
fit it, rather than continuous star formation. 

The columns are as follows: (a) Stellar Mass (1O 1O M0) (b) Stellar population age (Myr) from SED fitting, (c) Luminosity (xl0+ 28 ergs 
s — 1 Hz -1 ) from rest-frame UV continuum flux, (d) Star formation rate ( Mq yr - 1 ) estimated from rest-frame UV continuum flux, (e) 
Star formation rate ( Mq yr - 1 ) estimated from SED fitting, (f) Dust reddening derived from the SED fitting. 



3.4 Photometry, SED modelling and Stellar 
population properties 

We use broad-band photometry from the optical to the mid- 
infrared to determine the spectral energy distribution, and 
through fitting spectral evolutionary synthesis models we 
constrain the masses and ages of the stellar populations, 
and derive an independent estimate of the star rate to com- 
pare with that from the rest-frame optical nebular lines. We 
use the photometry from the VVDS catalog of the VVDS 
Deep Field 0226-0^] ( pVlcCracken et al.|2003||Le Fevre et aL 



2004) in the BVRI filters, obtained with the CFH12K cam- 



era on CFHT. For some galaxies this is supplemented by J 
and K s imaging with SOFI on the ESO NTT flovino et ah] 
2005 ; fTbmporin et al.||2008|). We supplement these magni- 



tudes with the recent CFHT Legacy Survey observations of 
this area obtained with MegaCam, using photomet ry in the 
griz from the TERAPIX CFHTLS-T0003 releas^l jSbert et 



al.|2006) . We note that all magnitudes used here are on the 



AB system. For the VVDS optical/near-IR photometry we 
are using MAG-AUTO (that is, "total" magnitudes derived 
by S Extract or). 

There is also Spitzer Space Telescope imaging of this 
field (the XMM LSS), obtained under the Spitzer Wide- 
area Infrared Extragalactic Survey (SWIRE, ISurace et al 



2005). We measured the IRAC photometry from version 4 



of the SWIRE team's reduced data product^] In the re- 
duced dataset, the IRAC pixels have been rebinned by a 
factor of 2 during the "drizzling" process from their original 
1//2 pix -1 , so the pixel scale is now 0'.'6. We convert the units 
of the mosaics from MJy/sr (the calibration from the Spitzer 
PBCD pipeline) to /iJy, requiring multiplying by 8.74. As 
with the optical/near-IR imaging, we work in AB magni- 



tudes, where 1 /iJy corresponds to 23.93 AB mags. Images 
in all four Spitzer channels we available (with central wave- 
lengths 3.6 /xm, 4.5 /im, 5.8 /xm & 8.0 /im), although the less 
sensitive channels 3 and 4 had no detections in the vicinity 
of our SINFONI observations, with 3 a limits of vtlab < 21.3 
and 21.1 at 5.8 /mi & 8.0 /xm respectively. For channels 1, 2, 
3 & 4 we used circular apertures of radius 2, 2, 2.5 & 3 pixels 
(2 / /4, 2 / /4, 3" h 3'/ 6 diameter) respectively, and an aperture 
correction of 0.7 mag (see |Eyles et al.||2005[ ). Once magni- 
tudes in each of the different wavebands had been obtained, 
the photometric data were then used to construct SEDs for 



each of our selected sources. We made use of the IBruzuall 
|fe Chariot] ( |2003| isochrone synthesis code (hereafter B&C), 
utilising the Padova-1994 evolutionary tracks (preferred by 
B&C). The models span a range of 221 age steps approx- 
imately logarithmically spaced, from 10 5 yr to 2 x 10 10 yr, 
although here we discount solutions older than ~2x 10 9 yr 
(the age of the Universe at z « 3). The B&C models have 
6900 wavelength steps, with high resolution (FWHM 3 A) 
and 1 A pixels over the wavelength range 3300 A to 9500 A 
and unevenly spaced outside this range. Throughout this 
paper, we adopt the |Chabrier| ( |2003| ) initial mass function 
(IMF). From the range of possible star formation histories 
(SFH) available, we focus on a constant star formation rate 
(SFR), as these Lyman-break galaxies are selected on their 
rest-UV continuum (i.e. have ongoing star formation). This 
star formation rate history is also used in the work of |Erb et] 



aLl( |2006b| ). For the constant SFR model, the B&C template 
normalization is an SFR of 1 Mq yr -1 . We also consider the 
possibility that the optical-infrared colours of objects within 
our sample could be due to intrinsic dust reddening, rather 
than an age-sensitive spectral break. We adopted the empiri- 
cal reddening model of Calzetti ( 1997 ), suitable for starburst 



galaxies. Table [4] summarise the properties of the stellar pop- 
ulation (absolute magnitudes, stellar mass, population age, 
1 Available from http://c encosw oamp.fr/VVDSphot /V VDS/ vvds.htnu] xtinction etc>) obtained for the ga l axies G f our sample. 
http://cencosw.oamp.tr/CEHTLS/ presents the CEHT-LS 



http: / / cencosw.oamp.fr/CEHTLS / presents 

magnitudes and photometric redshifts 

3 Available from http:/ /data. spitzer. caltech.edu/popular/swire/20061222_enhanced/XMM_LSS/irac/ 
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Figure 5. The spectral energy distribution for VVDS-3884, VVDS-8666, VVDS-7772 and VVDS-5183 from the broad-band CFHT 
optical/near-IR photometry, and Spitzer/IRAC channels 1 & 2. The best-fit continuous star formation model is overlaid. Bruzual & 
Chariot (2003) models have been used with a Chabrier (2003} IMF and solar metallicity. 

Table 5. Nebular emission properties 



Galaxy 


A neb (a 


) S gas (b) 


M gas (c) 




AGN? (e) 


12+log(0/H) (f) 


L(m ( g ) 


SFR , (h) 

neb v ' 


E(B - V) gas (i) 


SFR neb (j) 


VVDS-3884 


78.6 


521 


4.1 


0.76 


~ 40% chances 


8.57 ±0.02 


542 ± 27 


125 ±4 





125 


VVDS-8666 


70.2 


471 


3.3 


0.67 


no 


8.66 ±0.02 


429 ± 19 


97 ± 19 


0.35 


427 


VVDS-7772 


37.4 


<90 


<0.1 


<0.16 


yes 




< 29 (2a) 


< 5.1 (2a) 





< 5.1 


VVDS-5183 


74.6 


838 


6.3 


0.81 


~ 40% chances 


8.56 ±0.02 


949 ± 38 


231 ± 57 


0.40 


1257 


The columns 


are as 


follows: (a) 


area of the [O in] 


A5007 emission 


in kpc 2 , (b) gas surface density M@ yr 


_1 pc -2 , (c) Gas 


mass 



(1O 1O M0), (d) Gas fraction [i = M gas /(M gas + M*), (e) is the nebular spectrum dominated by an AGN ('no' implies it is dominated by 
star formation) (f) gas-phase oxygen abundance, (g) Luminosity (xlO 40 ergs -1 ), (h) star formation rate (MQyr -1 ) estimated from the 
H[3 luminosity (i) reddening suffered by gas, derived from SED fitting (j) dereddened star formation rate. 



3.5 Star Formation rates 

We computed star formation rates for the four galaxies in 
our sample deduced from the H/3, the SED fitting and from 
the rest-frame UV continuum emission. 

Ultraviolet-derived star formation rates were calculated 
from the broadband optical photometry, using /-band for 
VVDS-5183 at z = 3.7 and i?-band for the other galaxies at 
z ~ 3.3 to determine the UV continuum level around 1500 A. 
In the absence of dust, the UV continuum density per unit 
frequency (f u ) is approximately flat between the Lyman-a 
and the Balmer break for a constant star formation rate. 
The mean luminosity in frequency units is calculated using 
the following equation: 



Lisootergss" 1 Hz" 1 ) = 10 ( 



-OA-m(AB)) 



X363M0" 



We can then deduce SFRuv following Kennicutt (1998a) 
scaled to a |Chabrier| fl2003| ) IMF: 

SFR (M yr" 1 ) = 0.83 x 10" 28 L1500 (ergs s" 1 Hz" 1 ) (2) 

and the star formation rates derived from the rest-frame UV 
are shown in Table |4] 

We also derived star formation rates from the H/3 
emission-line luminosity. Given that no other Balmer line 
is observed, we are not able to compute the dust attenu- 
ation from the standard Balmer-decrement method. Thus, 
we decide to use instead the dust attenuation determined 



from the SED fitting. Followin g|Erb et al.|(|2006b|), we adopt 



E st ar(B — V) = E gas (B - V). Argence Lamareille (2009) 



have shown that when the dust attenuation is not computed 
from the Balmer-decrement method (e.g. with the observed 
Ha/H/3 ratio), the standard scaling-law which relates the 



emission-line luminosities and the star formation rates, i.e. 
SFR — ax L(line) ( Kennicuttj 1998a ), is not valid anymore. 
It has to be replaced instead by a power- la w relation, i.e. 



-(i) 



SFR — ax L(line) b with 6/1. We used thefCalzetti et al.| 



( 2000 ) extinction law. We thus combine equations 6 and 14 



of Argence & Lamareille ( 2009 ) to estimate the nebular star 



formation rates of our galaxies. We multiplied the results 
for star formation rate by a factor of 0.88 to convert from 
Kroupaj pOOl) IMF to |Chabrier| ( |2Q03| ) IMF. The resulting 
dust-corrected star formation rates from Hf3 are shown in 
Table [5] We found high values for SFR ne b corrected which 
should be taken with caution. It is indeed difficult to es- 
timate a reliable dust extinction from SED fitting at such 
high redshifts, where the relations between dust, age, and 
metallicity might be very different to what is observed and 
relatively well understood at low redshift. Thus, this values 
are indicative, and we use the non-dereddened values of the 
star formation rate, SFR^ ebl in the rest of the paper. 



3.6 Gas Masses and Gas Fractions 

We computed the gas mass using the empirical global 
Schmidt law, the correlation between star formation surface 
density and gas density ( 



Kennicutt 1998b). Assuming 



the galaxies in our sample follow such a law, we used the 
SFRn e b and the spatial extent of the [O ill] A5007 emission, 
an area A ne b (measured from the 2D intensity map of the 
[Oiii] A5007 emission, and deconvolved with the PSF, see 
|4.2| , to compute the star formation surface density (as is 
done in Erb et al. 2006b). From the star formation sur- 



face density, we deduced the gas density using equation 
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(4) in Kennicutt (1998b). We then deduced the gas masses 
by M gas = E gas A ne b. Using the stellar masses determined 
from the SED fitting, we can derive the gas fraction by 
\i — Mgas/ (Mgas + M*). Table [5] summarizes the results 
obtained for the four galaxies of our sample. All of these 
objects, except VVDS-7772, have high values of the gas frac- 
tion (/1 > 0.5), similar to that found in the z ~ 2 — 2.5 sample 
of |Law et ah] ( |2009l ) . 



4 PROPERTIES FROM SPATIALLY 
RESOLVED KINEMATIC 

4.1 Kinematic measurements 

We produced maps of the dynamics for the galaxies in 
our sample by using E3D, the Euro3D visualization tool 



(Sanchez et al. 2004), and code for the fitting and analy- 
sis of kinematics (e.g., Sanchez et al.| 2004 , 2005). Our main 
goal was to determine the kinematics of the ionized gas using 
the strongest emission line, [O ill] A5007 . For each spatial 
pixel (spaxel) we fit the [O ill] A5007 emission line region 
to a single Gaussian function, in order to characterize the 
emission line, and a pedestal to characterize any spectral 
continuum. We first smoothed the reduced cubes spatially 
with a Gaussian of FWHM = 3 pixels (0737). The line flux, 
FWHM, central wavelength and the continuum (pedestal) 
were then fitted. From the results of this fitting we obtained 
maps (see Fig. [6} of the [O in] A5007 emission line inten- 
sity, the relative radial velocity map (V r ) and the velocity 
dispersion (a). The dispersion map was corrected for the 
contribution of the instrumental dispersion, as determined 
from the FWHM of unblended and unresolved sky lines. Er- 
ror maps for the velocity and dispersion measurements were 
also computed, which are dominated by the effect of random 
noise in fitting the line profile, which produces large errors 
at low S/N. Errors ranged from 6 km s _1 to 24 km s _1 for 
the maps. 



4.2 Morphologic properties 

The left-most panel, (a), of Figure [6] shows the UV rest- 
frame morphology around 1500 A, with panel (b) showing 
the maps of the ionized gas morphology for the galaxies in 
our sample. For the four objects of the sample, the UV rest- 
frame morphology and the morphology of the ionized gas 
appear similar. VVDS-3884 presents a bright concentrated 
central region, with a diameter of 2.3 kpc. Two fainter re- 
gions slightly elongated toward the South and the Northwest 
surround the bright nucleus. VVDS-8666 exhibits an ellipti- 
cal morphology with a large peak in intensity at its centre in 
both the i?-band image and the line map. The galaxy VVDS- 
7772 looks compact, showing a peak in the ionized gas map 
slightly offset to the North-East. VVDS-5183 consists of a 
bright nucleus and a low surface brightness elongated region 
extended toward the North-East. 



For VV DS-8666 we used the GALFIT software ( |Peng 
et al.|200"2 ) on the CFHT R-b&nd image, to deduce the mor- 
phological parameters such as the center, the position an- 
gle, and the axial ratio. Assuming that the rest-UV light 
is dominated by stars in an inclined disk, we derived a po- 
sition angle of 6 = 41° East of North and an inclination 



of i = 51° (where i = 0° would be an edge-on disk). For 
the other galaxies in our sample, the images were too faint, 
compromised by a stellar diffraction spike, or lacked good 
nearby PSF stars for GALFIT to return satisfactory fits. 

We assume that the rest-UV morphology of the young 
stellar population will closely trace that of the nebular emis- 
sion from the ionized gas, and hence we infer the morpholog- 
ical parameters of the galaxies from the [O 111] map. We fit 
a 2D gaussian with the FWHM along the major axis, with 
the axial ratio and position angle as free parameters. A first 
guess of the position angle is estimated from the velocity 
map as the direction where the gradient in the velocity pro- 
file is maximum. It is then allowed to vary between ±15° 
in the rotation modelling (see |4.3[ ) which finally gives the 
best-fit position angle. If the galaxy is a disk, then the in- 
clination is estimated as cos _1 (6/a) where a is the radius of 
the major- axis (along the first guess of the position angle) 
and b the radius of the minor- axis (at 90°). We measured 
the area of the nebular emission (deconvolved with the PSF) 
from the two-dimensional [O ill] A5007 flux map for each of 
the galaxies (Table [5]). 



4.3 Rotation modelling 

Each galaxy is modelled by a pure, infinitely thin, rotating 
disk. The parameters of a model are the kinematic centre 
(xo,2/o), the position angle (PA), the inclination (z), the ve- 
locity offset (Vs) of the centre relative to the integrated spec- 
trum, and the velocity curve V c (r) where r is the radius from 
the kinematic centre. We also use the true physical veloc- 
ity dispersion ao as a model parameter, which we assume is 
constant, and represents the thickness of the rotating disk. 
The radial velocity V for any point is then defined with 
standard projection equations. Note that the position angle 
gives the direction of positive radial velocities. The veloc- 
ity offset accounts for redshift uncertainties. We define the 
velocity curve by two parameters, following Wri ght et al.| 
fl2007| ): the maximum velocity Vmax and the radius r c where 
this maximum velocity is achieved. In this model, below r c 
the velocity is defined as 



V c (r) = -xV m 
r c 



(3) 



The model computes V ro t, the asymptotic maximum rota- 
tion velocity at the plateau of the rotation curve, corrected 
for inclination (V ro t = Vmax/ sinz). 

In reality, the spatial resolution is limited by the seeing 
and the spaxel size. The observed radial velocity is thus the 
weighted convolution of the true radial velocity by the point 
spread function (PSF). The weights come from the flux map 
of the line used to compute the velocity map: a spaxel where 
the observed line flux is negligible will not contribute to the 
convolution. Additionally, the observed velocity dispersion 
accounts not only for the true physical dispersion, but also 
for the variations in the velocity field inside the width of the 
PSF. We generate a map of the velocity dispersion where 
the velocity gradient (determined from the modelling) has 
been subtracted. This yields a better measure of the intrinsic 
local velocity dispersion, rather than a raw value of a which 
may be inflated if the velocity gradient across a spaxel is 
large. From this map, we compute the flux- weighted mean 
velocity dispersion, cr mean (Table [6]). We also measure a ve- 



2D KINEMATICS AND PHYSICAL PROPERTIES OF z ~ 3 STAR-FORMING GALAXIES 9 




Figure 6. From left to right: (a) CFHTLS image (r-band for VVDS-3884, VVDS-8666 and VVDS-7772; i-band for VVDS-5183), (b) 
[O ill] A5007 flux map, (c) [O in] A5007 velocity field, and (d) [O in] A5007 velocity dispersion map obtained from Gaussian fits to the 
SINFONI data cubes after smoothing spatially with a two-dimensional Gaussian of FWHM = 3 pixels. The CFHTLS images and [O in] 
A5007 flux maps are color-coded with a linear scaling such that the values increase from light to dark. These data have been acquired 
with the 125x250mas sampling configuration of SINFONI. An angular scale bar of 3 /x , the redshift computed from the [Om] A5007 
emission line and the spatial rsolution are indicated for each galaxy. 



locity width for the emission lines from the extracted ID 
spectrum (i.e. the total integrated light); this am includes 
contributes both from the random motions and any ordered 
rotation or bulk motions. 

The beam smearing introduced by the PSF causes two 
significant effects. First, the observed velocity map appears 
smoothed, so the velocity gradient and the maximum veloc- 
ity are underestimated. Second, the observed dispersion map 
shows a peak of dispersion near to the kinematic centre. 

We can reproduce modelled velocity and dispersion 
maps by applying mathematically the same weighted convo- 
lution to the ideal velocity field. These modelled maps are 
compared to the observed maps by a \ 2 minimization. The 
parameters which minimize the \ 2 are computed in two suc- 
cessive grids with increasing resolution. First guesses for the 
kinematic centre, the position angle, and the inclination are 
set from the morphology. 



4.4 Properties of the ionized gas kinematics 

44.1 VVDS-8666 

The VVDS-8666 galaxy displays a well-resolved velocity gra- 
dient over no 5 kpc in projected distance with a peak-to-peak 
amplitude of 92kms -1 (uncorrected for inclination). There 
is some evidence of a flattening of the rotation curve, par- 
ticularly in the Eastern-most region of the galaxy. Taking 
into account the fact that these data are seeing limited, the 
smearing effect of the seeing decreases the peak-to-peak ob- 
served amplitude while mostly keeping the overall shape of 
the velocity shear. This shear appears to be aligned with the 
morphological major axis defined by the ionized gas, and is 
also nearly spatially coincident with the peak in the disper- 
sion map. This, combined with the fact that the nebular 
line flux distribution is centrally concentrated, leads us to 
conclude that the observed shear is consistent with a rotat- 
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Figure 7. Two-dimensional disk model after convolution with a 0.58 /x PSF for VVDS-8666 with the same spatial location as the observed 
velocity map (left). ID rotation curve (right; full diamond) extracted using an 'idealized' slit from the SINFONI [O ill] A5007 observed 
velocity map (see map (c) in Fig. [6j, overlaid with the best-fit model (solid line). 

major axis offset (kpc) 




-0.4 -0.2 0.2 0.4 
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Figure 8. Maps and curve are the same as in Fig.|7|but for the galaxy VVDS-7772. 



ing gaseous disk rather than a merger. Figure [7] shows the 
[O in] A5007 velocity map recovered from the rotation mod- 
elling (see Figure [6] panel c) and convolved with the PSF. It 
show also the one-dimensional relative velocity curves along 
the kinematic major axis for both the SINFONI observed 
velocity shear and the best-fit model. The velocity is well 
matched by our simple rotating disk model, leading to an 
inclination-corrected v ro t ~ 91=L26 kins" 1 . However, the ra- 
tio of the rotation velocity to the velocity dispersion is large 
(vrot/<Jmean = 1.52) compared with lower-redshift samples 
of galactic disks (v ro t/ crmean — 10 — 20, Dib et al. 2006) 



which probably indicates that this object is not mainly sup- 
ported by rotation, and there are significant random mo- 
tions. The presence of random motion is in particular con- 
firmed by a peak of ~ 50 km/s, located to the north of the 
galaxy, in the residuals of the velocity maps after subtraction 
of the rotating model. 



44.2 VVDS-7772 

The velocity map of VVDS-7772 (the companion galaxy 
of VVDS-8666) shows a smoothly varying gradient along 
the North- South axis, but without evidence for a flattening. 



However, this galaxy is very compact and only marginally 
resolved. We note that the velocity fields of both VVDS- 
7772 and VVDS-8666 are aligned along their displacement 
vector and close to the major axis of the resolved galaxy 
VVDS-8666, so this combined system might plausibly be 
part of a single rotating disk. The best fit simple rotation 
model give an asymptotic velocity of V ro t ^ 98 ± 45 km/s 
which is reached in ~ 5.4 kpc radius. The cr-map of the ion- 
ized gas of VVDS-7772 shows strong variation with posi- 
tion, as can be seen in panel (d) of Fig. [6] which prevents us 
from modelling the true physical dispersion; the weak con- 
straint is ao = =b 28kms _1 . In particular, a peaks around 
58kms _1 in the south-western edge of the galaxy, slightly 
offset with respect to the flux peak for which a decreases to 
around 50kms _1 , and drops further to ^30 km s _1 in the 
faint emission of the north-eastern region. However it is un- 
certain whether the peak of a ~ 58kms _1 represents a real 
feature or is simply noise in the extreme edges of the galaxy. 
Due to the compact size of this objects (r ~ 3kpc), it is dif- 
ficult to draw conclusions on the nature of the kinematics 
from our seeing-limited observations. 
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Figure 9. Maps and curve are the same as in Fig.|7|but for the galaxy VVDS-5183. 



44.3 VVDS-3884 

This galaxy presents an irregular velocity field, which is in- 
consistent with a smoothly varying velocity gradient along 
the major axis. It was therefore impossible to fit this dis- 
turbed shear with a single rotation model. The cr-map shows 
a peak in the local velocity dispersion displaced from the 
centre of the galaxy, again inconsistent with disk rotation. 
This galaxy has the largest velocity dispersion in our sample 
O -126kms _1 ). 



444 VVDS-5183 

This galaxy is dominated by an extended emission line com- 
ponent at the systemic redshift with a faint secondary ex- 
tension located ~ 2 kpc to the North-East. The two spatially 
distinct components are more apparent in the data cube (i. 
e. the velocity map) than in the [O in] A5007 intensity maps. 
Between the two components, the velocity seems to return 
to the systemic (central) velocity of the main component, 
i.e. the velocity of the whole system is not monotonically 
increasing with position. The faint secondary feature may 
represent a kinematically distinct star-forming region or a 
small galaxy in the process of merging with the brighter 
system. The primary component exhibits a smoothly vary- 
ing velocity field well fit by a rotating disk aligned with the 
morphological major axis. However the maximum of the ve- 
locity shear is small (V m ax ~ 30 km s _1 ), and once again this 
galaxy has significant random motion (high velocity disper- 
sion) with Vrot/cr mean — 0.95. 



4.5 Dynamical masses 

From the nebular line kinematics determined from our 3D 
data cubes, we are able to estimate the dynamical masses of 
the systems, subject to several caveats. Firstly, we assume 
that the nebular emission from the gas traces Keplerian mo- 
tions (i.e. the gas is not outflowing/ inflowing, and is in or- 
bit rather than unvirialized as might be the case during a 
merger) . Secondly, we cannot say with certainty that any of 
our galaxies are traditional rotationally-supported gaseous 
disks - the ground-based seeing precludes accurate measure- 
ment of morphological parameters (such as the crucial incli- 



nation angle). Thirdly, the measured velocity dispersion is 
significant when compared with any velocity gradient across 
these galaxies, implying that a putative disk may not be 
supported purely by rotation. 

We have measurements of the spatially varying velocity 
field (v) and velocity dispersion (a) from our 3D data cubes, 
and in the ID extracted spectrum we have a flux- weighted 
measure of the overall velocity spread in the galaxy (ctid) 
which arises from the combined effects of any systematic 
velocity shift across the galaxy and random motions. We 
can use this am as a crude estimate of the dynamical mass 
using the formula: 



Mdyn(criD) 



qas 

G 



(4) 



The factor C depends on the geometry of the system (in 
particularly the density profile): for a uniform sphere, C = 5, 
and Erb et al. (2006c) derive C — 3.4 for the more realistic 
scenario of a gas-rich disk with an average inclination angle. 
Using this value of C — 3.4, in Table [6] we present dynamical 
masses inferred from the g\d and the spatial extent of the 
[O in] A5007 line emission in our data cubes (r gas ) measured 
along the major axis and deconvolved with the seeing. 

Of course, the availability of 3D data cubes means we 
can model any spatially-resolved systemic velocity shift, as 



might arise in a rotating disk (see Section 4.3). For objects 
where the velocity shear is well fitted by the simple rotation 
model, we compute the dynamical masses using the formula: 



Md yn (Vrot) — 



V ro i T m odel 

G 



(5) 



where V ro t (the turnover of the rotation curve) has been 
inclination-corrected (a factor of sini, fitted by the model 
from an initial estimate based on the observed ellipticity ) . 
Both this asymptotic velocity, and the radius of the turnover 
(r model) are inferred from our model fits to the observed ve- 
locity maps to correct for the the effect of beam-smearing 
(due to the seeing). Only in one case (VVDS-8666) do we 
convincingly see this flattening of the rotation curve within 
the radius probed by our detected [O in] A5007 emission 
(Fig. [7), from which we infer M dyn = 2.9 x 10 10 M o . For 
galaxies VVDS-7772 and the primary component of VVDS- 
5183, the inferred turnover from the models lies beyond the 
range of data, and hence must be treated as unreliable (see 
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Table 6. Kinematics properties 



ID 


i(l) 


PA(1) 


Vmean (2) 


<?ID (3) 


V shear (4) 


Vrot (5) 


rgas 


Tmodel (6) 


M dyn (Vrot)/10 10 M Q (7) 


M dyn (a 1D )/10 10 M Q (8) 


VVDS-3884 


90±40 




72±29 


126±9 


50±12 




5.27±0.6 






65.7 


VVDS-8666 


51±25 


41±25 


60±16 


79±6 


46±16 


91±26 


4.89±0.6 


0.4±3 


2.9 


23.8 


VVDS-7772 


35±18 


73±51 


0±28 


53±12 


33±17 


98±45 


3.00±0.9 


1.4±2 


12.1 (> 0.76*) 


6.5 


VVDS-5183 


50±3lt 


-10±30t 


78±llt 


95±3 


28±6 


74±42 


5.84±0.4 (3.58+) 


4.3±3+ 


21.6+ (> 0.65*) 


41.3 



t estimated from the primary component only. 
* where we do not observe a flattening of the rotation curve, a robust lower-limit on the M^ yri is given by V^ hear r 9 as I ^ 
The columns are as follows: (1) morphological parameters: inclination 'i' and position angle 'PA'. (2) Flux- weighted mean velocity 
dispersions deduced from the rotation modelling, (3) Velocity dispersion of integrated spectrum from the [O ill] A5007fwhm. (4) Shear 
velocity f s hear — §(^max — %i n ) measured on the [Olll] A5007 velocity map, (5) asymptotic velocity inferred from the rotation 
modelling (6) radius inferred from the rotation modelling for which ^ s hear — V ro t (7) The total dynamical mass inferred from the 
velocity dispersion of the integrated spectrum (8) The total dynamical mass inferred from the rotation modelling 




• VVDS-7772 
A VVDS-8666 
■ VVDS-5183 
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logio o 1D /km s" 1 



Figure 10. Combined velocity, So.5 2 = 0.5V ro t 2 + o- me an 2 versus 
integrated linewidth am- The diagonal lines are the same as in 
|Weiner et ah] ( |2006| ), i.e. a 1 : 1 line and the |Rix et al.| ( |1997[ ) 
a = 0.6Vc line (V c is the circular velocity). Line width and So. 5 are 
correlated; the 0.5 pre-factor makes the combined velocity width 
a better estimate of velocity dispersion, so that the correlation is 
tighter and the galaxies closer to the 1 : 1 line. Our galaxies show 
good agreement with this. 



Fig. [8] and Fig. [^respectively); we can place a lower limit 
on the dynamical mass from the measured maximum ve- 
locity and spatial extent of the nebular emission: Md yn > 
Vshear r gas / G, which corresponds to Mdyn > 7 x 10 9 M in 
both these cases. The galaxy VVDS-3884 does not demon- 
strate evidence of significant ordered rotation, and we were 
unable to fit a simple rotating disk model in this case. 



5 DISCUSSION AND CONCLUSIONS 

We present the results of a study of the kinematic prop- 
erties of the ionized gas in four star-forming galaxies at 
redshift z ~ 3, from the spatially-resolved spectra of the 
[O 111] & H/3 nebular emission lines. We found that the 
ionized gas in these objects has high velocity dispersions 
(crmean ~ 60 — 70kms _1 ). We have also found that all our 



galaxies (except VVDS-7772) contain large quantities of gas 
as compared to their stellar mass and their dynamical mass 
inferred from rotation. With their high SFRs, it is likely that 
these objects are undergoing episodes of strong star forma- 
tion. Due to our seeing-limited observations, it is difficult 
to classify these objects with confidence as either 'disks' or 
mergers. In the case VVDS-8666 and its companion (VVDS- 
7772), we classify them as a "heated disk system" since for 
both galaxies the velocity structure of the ionized gas ap- 
pears to be consistent with rotation. In fact it is unlikely that 
these objects have a dynamically cold rotating disk of ion- 
ized gas, due to the significant velocity dispersions. VVDS- 
8666 has the highest metallicity, probably indicating a later 
evolutionary stage. The galaxy VVDS-3884 possesses a par- 
ticularly high nebular star formation rate as compared to its 
SED star formation rate which is consistent with this object 
experiencing a major recent burst of formation. This burst 
might originate from a major merging event, consistent with 
its complex kinematics. Finally, VVDS-5183 is mostly con- 
sistent with a merger on to a rotating disk. It has a huge star 
formation rate, a large fraction of gas and a low metallicity. 



Due to the uncertainties on the inclination angle for 
the objects in our sample, it is not possible to investi- 
gate the traditional Tully-Fisher relation using V ro t- How- 
ever, we can give an insight into a key question: Can we 
use the line widths of integrated nebular emission to get a 
good estimation of the dynamical masses of galaxies show- 



ing low values of v/a ? Figure 10 plots the combined veloc- 



ity scale So. 5 against the one-dimensional line width a id 
for the three galaxies for which we found a consistency 
with the presence of rotation. The combined velocity scale 
£0.5 2 — 0-5 Vrot 2 + cfmean 2 is estimated from the asymptotic 
velocity (V ro t) and the flux- weighted mean velocity disper- 
sion, both inferred from the rotation modelling (see Figure 
10). It represents an estimator, corrected from the smear- 



ing effect of the seeing (see section 4.3), of rotation velocity 



together with the presence of random motions. So. 5 is also 



well known to be strongly correlated with am (Weiner et 
aL|2006||Kassin et al.|2007| ). We found a very good correla- 
tion between S . 5 and a 1D for the VVDS-8666, VVDS-7772 
and the primary component of VVDS-5183. This shows not 
only that the one-dimensional line width a id can be used 
as a kinematic measure when estimating dynamical masses, 
but that Vrot alone cannot suffice to trace the majority of 
the dynamical mass. In fact, using only V ro t to investigate 
Tully-Fisher relation will push galaxies with cr mean > V ro t 
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X Law et al. (2009) 
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Figure 11. Maximum line-of-sight velocity shear vs. stellar 
masses for our four galaxies. The dashed line denotes (log(cr mea n)) 
of the |Law et al.| ( |2009| galaxies sample. We found that all 
our galax ies fall be l ow t his line. The solid line indicates the fit 
to the |Law et al.| ( |2009| ) sample (log (v s hear/kms _1 )= 0.30 x 
\og(M* /Mq) - 1.39). The box on this figure denotes the region of 
'SINS' galaxy population (see text). 



to erroneously low dynamical masses (see also Weiner et al. 



2006). 



Despite the small size of our sample, we can try to 
compare our results with those found by the 'SINS' survey 



( Forster Schreiber et al |2006| |Genzel et aI.|2006||Bouche et 



al.|2007||Shapiro et al.|2008| JGenzel et al.||2008| ), also using 
SINFONI on VLT, and the |Law et al.| fl2009| ) study using 
OSIRIS on Keck. Our four galaxies have a mean assem- 
bled stellar mass of (log (M*/M )) = 10.4, while the 'SINS' 
galaxy population has a larger mass of (log (M*/M©)) = 
10.9—11.0. More comparable to our sample are the 13 galax- 
ies detected by |Law et al. ( 2009| ), selected by the BX tech- 
nique, with stellar masses of (log(M*/M©)) « 10.1. We also 
found, in agreement with |Law et al.| ( |2009| ), that such galax- 
ies of intermediate mass at z ~ 3 possess extremely large 
quantities of gas in comparison to their stellar mass. The 
typical v/a for our sample is around ~ 0.4 — 1.5. The |Law| 
et al.| ( |2009| sample have ^0.8 whil e the 'SINS' gala xies 
have a range of ~ 2 — 4. In contrast to|Law et al. (2009), we 



have found evidence of AGNs in our sample. The connection 
between AGN and star formation could play an important 
role in shaping galaxies at z ~ 3 — 4. 

Figure [TT] plots the maximum observed shear velocity of 
galaxies, v s hear, as a function of stellar mass. Our z = 3 — 3.7 
objects appear to have values of ^ s hear similar to the high- 
est values detected by Law et al. ( |2009| but corresponding 
to the lowest values displayed by the galaxies of the 'SINS' 
survey. In fact at similar stellar masses our galaxies tend to 
have higher values of v s hear than Law et ah] (|2009| sample 
but similar values of v s hear to the Forster Schrei ber et al.| 
(2006) sample. We also found in general that the highest 



mass galaxies do not tend to have the highest velocity shear 
among our objects, similar to the Law et al. (2009} sample. 
Therefore, our sample at z = 3 — 3.7 seems to comprise in- 
termediate mass and intermediate shear velocity objects, in 



contrast to the sub-sample of Law et al. ( 2009| which seems 
to have negligible velocity shear and the more massive ob- 
jects, and the 'SINS' survey sample which have appreciable 
old populations where stable rotation dominates. 

The similarity in properties of our sample of galaxies at 
z ~ 3 - 3.7 with that of |Law et al.l (120091) at z ~ 2 - 2.5, 



seems to confirm that these high-redshift objects have not 
yet turned a significant fraction of their gas into a sizeable 
stellar population. Such galaxies also tend to possess high 
local velocity dispersions due to random motions intrinsic 
to the gas. If any velocity shear is present, it is not likely 
to be significant compared with the random motions, and 
probably occurs in a non- stable disk system. These galaxies 
also exhibit strong episodes of star formation and already 
possess metal rich gas. 

Drawing general conclusions about the nature and the 
properties of the dynamical structure of distant galaxies 
based on the current set of observations is obviously very 
challenging, given the small size of current samples and the 
spatial resolution limitations. However, initial insights into 
the kinematic properties of z ~ 2 — 4 galaxies on scales of a 
few kiloparsecs have been revealed by the recent IFU studies, 
although larger statistical samples and unbiased selection 
criteria are needed to give a quantitative picture of the mech- 
anisms at play in such objects. It already appears that the 
dynamical state of galaxies during this early period, where 
the star formation appears to be intense, cannot be easily 
classified based on what is observed at low redshift. Specif- 
ically, the star-forming galaxies at z ~ 3 (even those with 
measurable velocity gradients) exhibit large velocity disper- 
sions quite unlike the cold rotationally-supported disks seen 
at lower redshift. 



ACKNOWLEDGMENTS 

We would like to thank Aprajita Verma and Matthias Tecza 
for helpful discussions. We are very grateful to the VLT Ob- 
servatory for accepting this programme. The authors also 
thank Markus Hartung for help obtaining the observations 
during the observing runs and Sebastian Sanchez for pro- 
viding the code to create the kinematics maps from the 3D 
cube. The anonymous referee is greatly acknowledged for 
providing useful and constructive comments. The authors 
wish to recognize and acknowledge the significant contri- 
bution of the VVDS collaboration in providing the targets. 
In particular we thank Thierry Contini who helped greatly 
with the target selection. Part of this work was supported by 
the Marie Curie Research Training Network EuroSD; con- 
tract No. HPRN-CT-2002-00305. M. Lemoine-Busserolle is 
supported by the Gemini Observatory, which is operated 
by the Association of Universities for Research in Astron- 
omy, Inc., on behalf of the international Gemini partnership 
of Argentina, Australia, Brazil, Canada, Chile, the United 
Kingdom, and the United States of America. 



14 M. Lemoine-Busserolle et al. 



REFERENCES 

Argence B., Lamareille F., 2009, A&A, 495, 759 
Bouche, N., et al. 2007, ApJ, 671, 303 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Bunker, A. J., Moustakas, L. A., & Davis, M. 2000, ApJ, 
531, 95 

Bunker, A., et al. 2004, AN, 325, 139 
Calzetti, D. 1997, AJ, 113, 162 

Calzetti, D., Armus, L., Bohlin, R. C, Kinney, A. L., 
Koornneef, J., & Storchi-Bergmann, T. 2000, ApJ, 533, 
682 

Chabrier, G. 2003, PASP, 115, 763 

Davies, R. I. 2007, MNRAS, 375, 1099 

Dib, S., Bell, E., & Burkert, A. 2006, ApJ, 638, 797 

Eisenhauer, F., et al. 2003, SPIE, 4841, 1548 

Erb, D. K., Shapley, A. E., Steidel, C. C., Pettini, M., Adel- 

berger, K. L., Hunt, M. P., Moorwood, A. F. M., & Cuby, 

J.-G. 2003, ApJ, 591, 101 
Erb, D. K., Steidel, C. C., Shapley, A. E., Pettini, M., & 

Adelberger, K. L. 2004, ApJ, 612, 122 
Erb, D. K., Shapley, A. E., Pettini, M., Steidel, C. C., 

Reddy, N. A., & Adelberger, K. L. 2006a, ApJ, 644, 813 
Erb, D. K., Steidel, C. C., Shapley, A. E., Pettini, M., 

Reddy, N. A., & Adelberger, K. L. 2006b, ApJ, 647, 128 
Erb, D. K., Steidel, C. C., Shapley, A. E., Pettini, M., 

Reddy, N. A., & Adelberger, K. L. 2006, ApJ, 646, 107 
Eyles, L. P., Bunker, A. J., Stanway, E. R., Lacy, M., Ellis, 

R. S., & Doherty, M. 2005, MNRAS, 364, 443 
Faber, S. M., et al. 2007, ApJ, 665, 265 
Forster Schreiber, N. M., et al. 2006, APJ, 645, 1062 
Genel, S., et al. 2008, ApJ, 688, 789 
Genzel, R., et al. 2006, Nature, 442, 786 
Genzel, R., et al. 2008, ApJ, 687, 59 
Giavalisco, M., et al. 2004, ApJ, 600, L93 
Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 
Ilbert, O., et al. 2006, A&A, 457, 841 
Iovino, A., et al. 2005, A&A, 442, 423 
Kassin, S. A., et al. 2007, ApJ, 660, L35 
Kennicutt, R. C., Jr. 1998a, ARAA, 36, 189 
Kennicutt, R. C., Jr. 1998b, ApJ, 498, 541 
Kewley, L. J., & Elision, S. L. 2008, ApJ, 681, 1183 
Kroupa, P. 2001, Dynamics of Star Clusters and the Milky 

Way, 228, 187 

Law, D. R., Steidel, C. C, Erb, D. K., Larkin, J. E., Pettini, 

M., Shapley, A. E., & Wright, S. A. 2007a, ApJ, 669, 929 
Law, D. R., Steidel, C. C, Erb, D. K., Larkin, J. E., Pettini, 

M., Shapley, A. E., & Wright, S. A. 2009, ApJ, 697, 2057 
Lamareille F., et al., 2009, A&A, 495, 53 
Liang, Y. C, Yin, S. Y., Hammer, F., Deng, L. C, Flores, 

H., & Zhang, B. 2006, A&A, 652, 257 
Le Fevre, O., et al. 2004, A&A, 417, 839 
Le Fevre, O., et al. 2005, A&A, 439, 845 
Lemoine-Busserolle, M., Contini, T., Pello, R., Le Borgne, 

J.-F., Kneib, J.-P., & Lidman, C. 2003, A&A, 397, 839 
McCracken, H. J., et al. 2003, A&A, 410, 17 

Maiolino, R., et al. 2008, AfeA, 488, 463 

Modigliani, A., et al. 2007, |arXiv:astro-ph/0701297| 
Nesvadba, N. P. H., Lehnert, M. D., Davies, R. I., Verma, 

A., & Eisenhauer, F. 2008, A&A, 479, 67 
Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 
Peng, C. Y., Ho, L. C, Impey, C. D., & Rix, H.-W. 2002, 



AJ, 124, 266 

Rix, H.-W., Guhathakurta, P., Colless, M., & Ing, K. 1997, 

MNRAS, 285, 779 
Salpeter, E. E. 1955, ApJ, 121, 161 

Sanchez, S. F., Becker, T., & Kelz, A. 2004, Astronomische 

Nachrichten, 325, 171 
Sanchez, S. F., Garcia-Lorenzo, B., Mediavilla, E., 

Gonzalez- Serrano, J. I., & Christensen, L. 2004, ApJ, 615, 

156 

Sanchez, S. F., Becker, T., Garcia-Lorenzo, B., Benn, C. R., 

Christensen, L., Kelz, A., Jahnke, K., & Roth, M. M. 2005, 

A&A, 429, L21 
Shapiro, K. L., et al. 2008, ApJ, 682, 231 
Smith, J. K.., et al. 2004, MNRAS, 354, L19 
Stark, D. P., Swinbank, A. M., Ellis, R. S., Dye, S., Smail, 

I. R., & Richard, J. 2008, Nature, 455, 775 
Surace, J. A., Shupe, D. L., Fang, F., Evans, T., Alexov, 

A., Frayer, D., Lonsdale, C. J., & SWIRE Team 2005, 

Bulletin of the American Astronomical Society, 37, 1246 
Temporin, S., et al. 2008, A&A, 482, 81 
Tremonti, C. A., et al. 2004, ApJ, 613, 898 
van Starkenburg, L., van der Werf, P. P., Franx, M., Labbe, 

I., Rudnick, G., & Wuyts, S. 2008, A&A, 488, 99 
Weiner, B. J., et al. 2006, ApJ, 653, 1027 
Williams, R. E., et al. 1996, A J, 112, 1335 
Wright, S. A., et al. 2007, ApJ, 658, 78 
Wright, S. A., Larkin, J. E., Law, D. R., Steidel, C. C, 

Shapley, A. E., Erb, D. K. 2009, ApJ, 699, 421 



